Vortex Penetration In Magneto-Superconducting Heterostructures. 
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We report our results on vortex penetration in two realizations of heterogeneos magneto- 
superconducting systems (HMSS) based on London approach; semi-infinite ferromagnetic(FM)- 
superconducting(SC) bilayers and a FM dot on a semi-infinite SC film. In the first case, we study 
quantitatively the vortex entry in FM-SC bilayers which manifests Bean-Livingston-like vortex bar- 
rier, controlled by FM film's magnetization m and SC film's Ginzburg parameter k. In the second 
case, we investigate the conditions for sponteneous vortex creation and determine the position of 
vortex for various values of magnetization and the dot's position. 
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PACS Number(s): 74.25.Dw, 74.25.Ha, 74.25. Qt, 74.78, 



Heterogeneous magneto-superconducting systems 
(HMSS) are made of ferromagnetic (FM) and supercon- 
ducting (SC) pieces separated by thin layers of insulating 
oxides. In contrast to the case of a homogeneous ferro- 
magnetic superconductor studied during the last two 
decades, the two order parameters, the magnetization 
and the SC electron density do not suppress each other 
[1,2]. In HMSS, the strong interaction between FM and 
SC components stems from the magnetic fields generated 
by the inhomogeneous magnetization and the supercur- 
rents as well as SC vortices. Strong interaction of the 
FM and SC systems not only gives rise to a new class of 
novel phenomena and physical effects, but also shows the 
important technological promise of devices whose trans- 
port properties can be easily tuned by comparatively 
weak magnetic fields. 

Various theoretical realizations of HMSS have been 
proposed by different groups, such as arrays of mag- 
netic dots on the top of a SC film [1,3], ferromag- 
netic/superconducting bilayers (FSB) [4], and magnetic 
nanorods embedded into a superconductor [5], whereas 
only sub-micron magnetic dots covered by thin SC films 
have been prepared and studied [6-8]. The experimen- 
tal samples of FM-SC hybrid systems were prepared by 
means of electron beam lithography and lift-off tech- 
niques [9] . Both in-plane and out-of-plane magnetization 
was experimentally studied. The dots with magnetiza- 
tion parallel to the plane were fabricated from Co, Ni, 
Fe, Gd-Co and Sm-Co alloys. For the dots with magne- 
tization perpendicular to the plane which requires high 
anisotropy along hard-axis, Co/Pt multilayers were used 
[10]. 

In the most of theoretical studies, SC subsystem is 
considered to be infinite size for the sake of computa- 
tional simplicity. To this date, there has not been an 
analytical analysis of boundary and edge effects in FM- 
SC heterostructures, though vortex entry conditions in 



type II superconductors are previously studied [11-13]. 
However, from both experimental and theoretical point 
of view, finite or semi-infinite systems are more inter- 
esting and realistic, and their study offers better under- 
standing of vortex matter in HMSS. Author also believes 
that analytical and quantitative study of aforementioned 
systems will shed light on solving other open problems 
pertaining to HMSS. For example, we earlier predicted 
that in a finite temperature interval below the SC tran- 
sition the FSB is unstable with respect to SC vortex 
formation in FM-SC bilayers (FSB) [14]. The slow de- 
cay (oc 1/r) of the long-range interactions between Pearl 
vortices makes the structure that consists of alternating 
domains with opposite magnetization and vorticity ener- 
getically favorable. It is possible that the long domain 
nucleation time can interfere with the observation of de- 
scribed textures. We also expect that domain nucleation 
starts near the edge, which makes qualitative study of 
edges in aforementioned systems necessary. Quantitative 
study of this dynamic process is still in progress. For this 
purpose, and having been motivated by current interest 
in HMSS, in this work, we attempt to study vortex en- 
try conditions in HMSS. To our purpose, we work with 
a method based on London-Maxwell equations, which is 
fully explained elsewhere [15]. London approach works 
well for large Ginzburg (k = A e ///£ » 1 ) parameter, 
where A e // = \\/d is the effective penetration depth 
[16], and £ is coherence length. Indeed, for thin SC films, 
Ginzburg parameter is on order of 50-100. Previously, 
our method was introduced for vortex structures in in- 
finite films. Here, we extend it to semi-infinite systems. 
To this end, we benefit from Kogan's work on a Pearl 
vortex near the edge of SC thin film in which SC piece's 
size is considered to be semi-infinite [17]. Likewise, we 
consider FM subsystems on semi-infinite SC and FM sub- 
systems in which, we assume that magnetization points 
perpendicular to the FM film's plane. 
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In this work, wc first consider semi-infinite SC and 
FM films and study vortex entry barrier. Our calcula- 
tions show that there exists Bean-Livingston-like surface 
barrier [18] for the vortices created by FM film. Next, we 
consider a circular magnetic dot near the film's edge and 
investigate the conditions for vortices to appear and their 
configurations. It turns out that, in contrast to the infi- 
nite systems, vortices are not trapped right at the dots 
center, but they are shifted slightly from the center to the 
SC film's edge or opposite direction, depending on the 
dot's magnetization, position and size. Physics behind 
this effect is simple. In the semi-infinite systems, vortex 
interacts with both its image vortex and the magnetic 
dot. The competition between these two attractions de- 
termines vortex's position. The outline of this articles is 
as follows: in the first section, we introduce the method 
to study edge effects in FM-SC systems. In the next 
section, we apply our method based on Maxwell-London 
equations to two different cases; semi-infinite FM-SC bi- 
layers and FM dot on a semi-infinite SC film. We con- 
clude with results and discussions. 



We start with London equation for the vortices with vor- 
ticity rij located at r^, 

h+-^V xj, = c/) zJ2 n j 6 ( r (1) 



where j s is supercurrent density in the SC subsystem. In 
the presence of FM subsystem which is also considered to 
be very thin and located at x-y half plane as SC subpiece, 
Eq.(l) turns to 

h + x j = x Jm + <p z 2^ n 3 8{v - r 3 ), 



(2) 



where j = j s + j ro and j m = cV x m. Averaging Eq.(2) 
over the thickness of SC film, one finds 



I. METHOD 



(3) 



Finite and semi-infinite systems are not as easy and 
straightforward as infinite systems, and they usually re- 
quire more careful treatment due to the boundary of the 
systems. Earlier, Kogan developed a clever technique 
based on London approach to study a vortex near the 2d 
film's edge [17]. While developing our method, we stick 
to his technique and geometry in which a very thin SC 
film is located at x-y half plane while its edge is at x=0 
(see Fig.l), and generalize his method for more than one 
vortex. We also assume that no vortex is closer to the 
SC film's edge than coherence length £, because London 
theory does fail in the vicinity of £. 
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FIG. 1. Semi-infinite FM-SC bilayer. 



where g is the 2-d current density, which can be cal- 
culated by solving Eq.(3) together with the Biot-Savart 
integral equation and the continuity equation V • g = 0. 
In terms of the surface current, the Biot-Savart equation 
is given by 



h z = -f dV[g(r')xR/i? 3 ], 



(4) 



where R = r — r' . Defining the surface current density in 
terms of a scalar function G(r) as g = V x G(r)z, using 
R./R 3 = V(l/R) and integrating Eq.(4) by parts, one 
can obtain 



h z = - [ < 

c Jx'>0 



R 



dy' 



, ( d x >G{r') 
R 



x'=0 



(5) 



where the first term gives the contribution from the en- 
tire surface current distribution whereas the second term 
is the contribution from the film's edge. The direct sub- 
stitution of Eq.(5) into Eq.(2) gives 



lx'>0 
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^G(r') 
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+ 47rAV 2 G(r) = -c^ ^ n A v - r i) - 4ttA(V x g m ) 2 
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(6) 



Solving the above equation in half plane is difficult. How- 
ever, this difficulty can be removed by solving (6) in the 
Fourier space and using the boundary conditions. That 
is, at the Alms's edge ( x = ), the normal component 



of current density is zero, namely g x (0,y) = 0, whereas, 
at infinity the current distribution vanishes. This implies 
that the scalar function is constant at the films bound- 
aries. For simplicity, it can be set to zero. To have G 
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The asymptotics of the magnetic field at small and large where E v is the vortex energy, E vm is the interaction of 
distances are vortex and magnetic subsystem and finally E m is self- 

energy of magnetic subsystem, which will be ignored at 
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4A |r + a| 



r << A (18) 



further calculations, since it is inappropriate for our prob- 
lem. Vortex energy is calculated by Kogan [17] as 



r >> A. 



(19) 



|r + a| 3 

The total energy of FM-SC system reads 
E E v -\- E vvtl -\- E m , 

8A 7T / Xi \ t—^ 



\ - n t <p 



(21) 



where ^ denotes the restricted sum in which only i = j 
(20) and i > j are taken into account. Eq.(21) leads to 
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(22) 



where &o(x) — Yq(x) — H^(x). Vortex- magnetization in- 
teraction energy is calculated as in [15]: 

= J • a m d 2 x -\f™- b v d 2 x, (23) 

where integration is performed over the half space. Note 
that we take A/£ = 50 in our numerical calculations. 



A. Semi-Infinite FM-SC Bilayers 

In this part, we study a semi-infinite FM film on top of 
a semi-infinite SC film. Both films are taken to be very 
thin, lie on x-y half plane whereas their edges are located 
at x = (see Fig.l). We assume that FM film has uniform 
magnetization along z direction and has high anisotropy, 
so that its magnetization does not change direction due 
to magnetic field of vortex. The magnetization of FM 
film reads 



m = m0(x)8(z)z. 



(24) 



Magnetic current in real space and Fourier space is given 
as 



-mcS(x)y, 



-2irmc6(ky)k y . (25) 



Substituting Eq.(24) into Eq.(10) , one can find the scalar 
potential as 



G m (k) = SirXmci 



k x 5(ky) 



(26) 



fc(l + 4AA;)' 

Taking the inverse Fourier transform of Eq.(26), we find 

where f(x) = dk x sm(k x x)/(l + k x )). The asymp- 
totics of f(x) are 



+ x(\n(x) + C - l),x « 1, 
x » 1. 



(28) 



Using Eq.(12) and Eq.(13), z component of the screened 
magnetic field at z = due to FM film, reads 



k x sin(fca;^-) 
1 + 1**1 X ' 



(29) 



The magnetic field decays as 1/x for x « A and 1/x 2 for 
x » A. In order to study vortex configuration, we need 
to calculate total effective energy of the system. To this 
end, we consider a simple case, namely a vortex with a 
single flux located at r = a. For this case, vortex energy 
for a single vortex reads (see Eq.(22)) 



E v = 



16tt 2 A 



8A 7T / a \ 

In — ; $ — 

e A f 2 V2A/ 



(30) 



wheras the vortex-magnetization interaction energy can 
be calculated by means of Eq.(23). 



E vm = -m<j>o 



!--/(- 

7r V4A/ 



(31) 
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FIG. 2. Phase diagram of a single vortex created by 
semi-infinite FM film. In the region below curve, vortex does 
not appear, whereas it becomes energetically favorable in area 
above the curve. 



The sum of Eq.(22) and Eq.(31) gives the effective to- 
tal energy of the system. Vortex becomes energetically 
favorable when effective total energy becomes less than 
zero. Equating the effective energy to zero, one can ob- 
tain the curve for spontenous creation of vortex. This 
curve, 




FIG. 3. The effective energy versus the vortex's position. 
When m < m c i, vortex does not appear. 



_H%)-l(H,{f- x )-Y {f- x )) 
1 - If (fx) 



scperates the regions where the vortex appears sponta- 
neously and does appear as seen in Fig. 2. For large 
values of to^o/co ratio, the vortex comes out near the 
edge. On the other hand, it prefers going further away 
from the surface for small ratio of m<fro/£o- We can es- 
timate the minimum value of magnetization of FM film 
through effective energy for infinite films [15], which is 
E e ff — £o ln(A/£) — m</>o- Equating this equation to zero 
and solving it for m, we find m c i = </>o/(167r 2 Aln(A/£)). 
When magnetization exceeds this value, the vortex ap- 
pears very far away from the edge. In order to get vor- 
tex appear close to the edge, m must be significantly 
larger than m c \. Another interesting thing is that the 
system manifests Bean-like surface barrier for the vor- 
tex. The surface barrier is controlled by m</> /eo and 
Ginzburg parameter n. We analyze three regimes for 
this ratio for fixed k = A/£. When m < m c i, vor- 
tex does not appear (see Fig. 3). In the second regime, 
m c i < m < m C 2, vortex prefers going further away from 
the surface, whereas, when m > to c2 , the barrier dis- 
appears (see Fig. 4). m c2 is the second critical magne- 
tization, at which barrier disappears, and can be calcu- 
lated through the condition that the slope \dE to t/d%\x=£ 
is zero, which gives m0 o /eo ~ 27tk/1ii(4k). When ra- 
tio rrufio/eo is greater than this , the barrier disappears. 
Physically, two contributions play an important role for 
the vortex barrier. Namely, the vortex is attracted to SC 
film's edge through its attraction towards image vortex 
whereas it is repelled by FM film's edge. Competition 
between these two factors controls the barrier. 
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FIG. 4. The effective energy versus the vortex's position. 
When m c i < m < m C 2, the surface barrier shrinks toward the 
edge of SC film, and vortex is created little further from the 
edge. 
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FIG. 5. The effective energy versus the vortex's position. 
When m > m C 2 barrier disappears, and vortex can be seen 
anywhere in SC film. 
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B. FM Dot on a Semi-Infinite SC Film 



FM disc located at r<j = (xd, 0), 



In this case, we study a circular FM disc on top of 
semi-infinite SC film. Earlier, we studied the conditions 
for the vortex states to appear on a similar system, in 
which SC film, however was infinite. Due to the circular 
symmetry of FM dot, vortex was appearing at the dot's 
center. In this section, we study the vortex states in more 
realastic case and investigate the role of edge effects on 
the spontaneous formation of vortex due to FM dot. To 
this end, we start with the magnetization of a circular 



m = m8(R — |r — rd\)S(z)z, 



(33) 



where R is the radius of circular dot. The Fourier trans- 
form of Eq.(33) reads 



nik = 2nmR 



J 1 {kR)e 



-ik x x d 



k 



(34) 



Using Eq. (10,12,13) together with Eq.(34), one can cal- 
culate the screened magnetic field due to FM dot as 



h z (r, z) — AnmXR 



Ji(fcfl)[Jo(fc|r - r d |) - J (fc|r + r d \)]k 2 e 
1 + 4Afc 



-k\z\ 



-dk. 



(35) 



Outside the dot, magnetic field decays rapidly in space, From Eq.(23), the vortex-magnetic disc interaction en- 
namely, for r « A, ~ 1/r 3 , whereas for r « A, <~ 1/r 5 . ergy reads 

tp a u f°° t (un\ lM k \ r d - a|) - J (k\r d + a|)l 

E vm = -m(p Q R Ji{kR) -dk. (36) 

Jo 1 + 4Afc 



After we formulate the total effective energy as E = 
E v + E vrn , where E v is given in Eq.(30), we study the 
conditions for a vortex to appear spontaneously. The 
criteria for spontaneous vortex formation is that effec- 
tive energy becomes negative. However, vortex in semi- 
infinite systems also interacts with its image. Therefore, 
it is necessary to minimize total effective energy with re- 
spect to the vortex position. To this end, we first fix the 
dot's location and value of to</>o/£o and vary the vortex's 
position afterwards, to find the minimum total effective 
energy. In our calculations, we investigate where vortex 
first comes out, and how it is shifted with further increase 
of rn</>o/eo- For this purpose, we determine vortex's po- 
sition for different values of the r<j/A, R/X and m<j) /eo, 
by optimizing the total effective energy. Our results are 
shown in Tabled and Fig. 6. 

According to our calculations, vortex appears first 
close to the edge except r^/R > 10 case, in which it sits 
at the dot's center. On the other hand, when m^o/eo 
is increased further, vortex is first shifted towards the 
dot's center. With further increase of mcpo/so, it drifts 
away from the dot's center. However, this is not always 
general picture. In the case of rj/R > 5/3, vortex is lo- 
cated at the dot's center even for large values of m<fo/eo- 
However, for larger dot's sizes (ra/R ~ 1 and r^/A > 2) 



vortex first appears away from the dot's center. This sit- 
uation differs from the vortex in infinite system. In the 
latter, only force acting on vortex stems from the vortex- 
magnetization interaction, and due to the dot's circular 
symmetry, it comes out at the dot's center. However, in 
this case, there is another force coming from vortex-image 
vortex, which decays slow ~ 1/r for large distances r > A 
and pulls vortex towards the SC film's edge, whereas the 
force exerted by the magnetic dot pushes vortex towards 
the dot's center. As a result, vortex's position is deter- 
mined by the balance between these two forces. 



TABLE I. The position of vortices for different values of 
the 7\j/A, R/X and m0o/eo- The two columns on the left are 
input. 
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R/X 
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2.0 


2.04 
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2.0 


2.32 


187 


3.0 
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3.12 
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3.0 


3.44 


132 
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4.0 


4.20 


11 


4.0 


4.0 


4.56 


185 
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R/A-0.3 
R/A-0.5 
R/X=0.8 



0.42 0.44 0.46 0.48 0.5 0.52 0.54 0.56 0.58 0.6 

all 



0.85 

all 



(a) r d /X = 0.5 



(b) r d /X = o.i 
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FIG. 6. Position of vortex a/A versus m^o/eo for various locations r d j\ and sizes R/\ of magnetic dot 



II. CONCLUSIONS 

In this article, we studied vortex entry conditiions in 
HMSS. First, we generalized Kogan's method for quan- 
titative study of semi-infinite HMSS. For applications, 
we first considered semi-infinite FM film on top of semi- 
infinite SC film. The quantative analysis of this system 
showed that the vortex undergoes Bean-like barrier which 
is controlled by two intrinsic properties of system; FM 
film's magnetization m and Ginzburg parameter k. Note 
that our result is valid only for a single flux. Our study 
for the case of several fluxes will be published elsewhere. 
Secondly, we studied a single circular FM dot on a semi- 
infinite SC film. We analyzed the conditions for spon- 
teneous vortex creation and vortex location for various 
positions of the dot. It turns out that the vortex does 
not always appear at the dot's center, which differs from 
the case in which similar dot on an infinite SC film. In 
closing, there are two important contributions in semi- 
infinite HMSS; attraction of vortex to the edge through 
its image vortex and vortex-magnetization interaction. 
As a result of competion betweeen these two factors, pe- 



culiar physical effects which do not come out in infinite 
HMSS, appear. In this work, we studied the simplest 
cases to get idea about edge effects in HMSS. However, 
there are still several interesting realizations that can be 
studied via the method that is developed here. We leave 
them for possible future works. 
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